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DELAY- LOCK  TRACKING  OF  BINARY  SIGNALS 
J.  J.  Spilker,  Jr.  ,*  Member,  IRE 


Summary 


This  paper  presents  the  theory  of  operation  and  an  evaluation  of  performance  of  a  delay- 
lock  tracking  system  for  binary  signals.  The  delay-lock  discriminator  is  a  nonlinear 
feedback  system  which  employs  a  form  of  cross-correlation  in  the  feedback  loop.  Its 
function  is  to  estimate  continuously  the  relative  delay  between  a  reference  signal  and 
a  delayed  version  of  that  signal  which  is  perturbed  with  additive  noise.  Binary  maximal- 
length  shift-register  sequences  are  used  as  the  signal  because  they  can  easily  be  re¬ 
generated  with  any  desired  delay  and  because  they  possess  desirable  autocorrelation 
functions . 

Problems  of  target  search  and  acquisition  are  studied.  The  system  performance  in 
the  presence  of  additive  Gaussian  noise  is  discussed.  Computations  are  made  of  the 
effect  of  amplitude-limiting  the  received  data  on  the  system  noise  performance. 


♦Electronic  Sciences  Laboratory,  Lockheed  Missiles  &  Space  Company,  Palo  Alto 
California.  ’ 
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DELAY-LOCK  TRACKING  OF  BINARY  SIGNALS 
J.  J.  Spilker,  Jr.,  Member,  IRE 

I  INTRODUCTION 

1  Q 

The  delay-lock  discriminator  has  been  described  previously  ’  as  a  device  for  track¬ 
ing  the  delay  difference  between  two  correlated  waveforms.  The  discriminator  is  a 
nonlinear  feedback  system  which  employs  a  form  of  cross-correlation  in  the  feedback 
loop.  This  device  tracks  the  delay  of  a  broadband  signal  much  in  the  same  manner  that 
a  phase-lock  discriminator  tracks  the  phase  of  a  sinusoidal  signal.  In  this  paper,  a 
modified  version  of  the  delay-lock  discriminator  is  described  which  is  particularly 
designed  to  track  binary  signals  generated  from  feedback  shift-registers.  One  of  the 
main  reasons  for  the  interest  in  this  type  of  signal  is  that  it  is  easy  to  regenerate  the 
signal  in  the  discriminator  for  use  in  the  cross-correlation  operations  with  any  desired 
amount  of  delay.  Furthermore,  certain  classes  of  these  sequences  possess  desirable 
autocorrelation  properties;  e.g.  ,  the  maximal -length  linear  shift-register  sequences3,4 
have  "two-level"  autocorrelation  functions  and  can  be  designed  to  have  periods  which 
are  long  enough  so  that  ambiguities  are  of  no  concern. 


J.  J.  Spilker,  Jr.,  D.  T.  Magill,  "The  delay-lock  discriminator  -  an  optimum  track¬ 
ing  device,"  Proc.  IRE,  vol.  49,  pp.  1403-1416,  September,  1961. 

M.  R.  O'Sullivan,  "Tracking  systems  employing  the  delay -lock  discriminator,"  IRE 
g Trans.  ,  vol.  SET-8,  pp.  1-7,  March,  1962. 

B.  Elspas,  "The  theory  of  autonomous  linear  sequential  switching  networks,"  IRE 
^ Trans.  ,  vol.  CT-6,  pp.  45-60,  March,  1959. 

S.  W.  Golomb,  L.  R.  Welch,  "Nonlinear  shift-register  sequences,"  Jet  Propulsion 
Laboratory,  California  Institute  of  Technology,  Pasadena,  California,  Memorandum 
#20-149,  October,  1957. 
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In  an  actual  tracking  operation,  the  binary  sequence  is  applied  to  an  RF  carrier  and  is 
transmitted  to  the  target,  where  it  is  returned  to  the  tracking  equipment  either  as  a 
radar  reflection  or  via  a  transponder.  The  RF  signal  is  then  demodulated,  and  the 
demodulated  binary  sequence  plus  additive  noise  is  fed  into  the  delay-lock  discrimina¬ 
tor  to  obtain  the  estimate  of  the  delay  between  the  original  and  the  received  sequences. 
This  paper  is  concerned  only  with  the  delay  estimation  operation  after  demodulation  of 
the  RF  carrier. 

The  analysis  of  this  feedback  system  is  presented  in  four  main  parts: 

1.  Derivation  of  the  system  equations 

2.  Evaluation  of  noise  effects, caused  by  receiver  noise  and  system 
self-noise 

3.  Determination  of  transient  performance  —  acquisition  and  pull-out 
effects 

4.  Computation  of  noise  effects  with  amplitude-limiting  inputs 

II  DESCRIPTION  OF  THE  DISCRIMINATOR 

A  block  diagram  of  the  modified  delay -lock  discriminator  is  shown  in  Fig.  1.  The  re¬ 
ceived  signal  originated  from  an  n-stage  maximal -length  linear  feedback  shift-register 
(FSR),  and  has  a  period  MA  =  (2n  -  1)A  ,  where  A  is  the  digit  width.  This  binary 
signal  plus  additive  Gaussian  white  noise  is  fed  into  the  cross -correlation  network, 
where  a  comparison  is  made  with  time-displaced  versions  of  the  same  pseudo-random 
binary  signal  as  used  in  the  transmitter.  As  is  shown  in  a  later  section,  the  implemen¬ 
tation  of  this  cross -correlator  can  be  substantially  simplified  if  the  received  data  is 
first  converted  to  binary  form  by  means  of  an  amplitude-limiter.  Then  the  multipli¬ 
cation  of  ±1  binary  digits  can  be  replaced  by  a  modulo  2  adder. 
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The  output  of  the  cross-correlator  network  in  Fig.  1  is5 

k<5s  [t  +  T(t)]  j/Ts[t  +  T(t)]  +  n(t)j  (1) 

where  6s  (t)  =  s(t  +  A)  -  s  ( t  -  A)  ,  and  Pg  is  the  received  signal  power.  As  will 
be  shown,  this  product  contains  a  low-pass  spectral  component  which  serves  to  keep 
the  discriminator  accurately  tracking  the  target  delay  once  the  system  has  been  "locked 
oa"  The  low-pass  filter  which  follows  the  cross-correlation  network  is  designed  on 
the  basis  of  the  expected  dynamics  of  target  motion,  and  its  purpose  is  to  remove  as 
much  noise  and  other  interference  as  possible.  The  output  of  the  loop  filter  is  an  esti¬ 
mate  of  the  time  derivative  of  delay  (proportional  to  the  target's  radial  velocity)  and  is 
used  to  control  the  clock  rate  of  the  FSR. 

The  delay  estimate,  T  ,  can  be  easily  obtained  by  sensing  the  time  instants  when  the 
transmitter  and  receiver  FSRs  go  through  a  particular  state,  e.g.  ,  the  "all  ones" 
state,  and  by  determining  the  time  difference  between  these  time  instants. 


Ill  SYSTEM  EQUATIONS 

In  order  to  analyze  the  response  of  the  system,  the  signal  cross -correlation  term  in 
(1)  is  expressed  as 


k6s  (t  +  T)  s  (t  +  T 


/F  £ 

v  s 


°A<e 


ns(t,e) 


(2) 


A  A 

where  e  -  T  -  T  ,  is  the  delay  error.  The  term  which  is  not  explicitly  dependent 
on  time,  namely ,  D^(e)  ,  the  discriminator  characteristic,  has  been  separated  from 


g 

Notice  that  6s/2A  is  reminiscent  of  the  expression  for  the  differentiated  signal  which 
is  used  in  the  cross -correlation  operation  in  footnote  reference  1: 

ds  _  lim  s(t  +  A)  -  s(t  -  A) 

'  dt  A— 0  2  A 
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the  remainder,  n  (t,e)  ,  the  self-noise  term.  The  discriminator  characteristic  can 
®  6 

easily  be  shown  to  be 


DA(e)  =  E  [6s(t  +  T)  s  ( t  +  T)] 
-  Rg  ( c  -  A)  -  Rs(e  +  A) 


(3) 


where  R  (  )  is  the  signal  autocorrelation  function. 

In  this  paper  the  binary  signal  has  amplitudes  s  =  ±1.  The  autocorrelation  function7,8 
and  the  discriminator  characteristic  for  this  maximal -length  sequence  are  plotted  in 
Fig.  2a,  b.  As  can  be  seen,  0^(0  varies  linearly  with  e  for  |e|  <  A  ,  and  is 
zero  for  2 A  <  |ej  <  (M  -  2)  A  .  Since  s(t)  is  periodic  with  period  MA  ,  DA(e) 
has  the  same  periodicity. 


The  self-noise  term  in  the  output  of  the  multiplier  network  can  be  expressed  as 
ng(t,e)  =  s(t+A  +  f)s{t+T)-s(t-A+T)s(t+T)-DA(e) 

A  A  (' 

=  S(t  +  A  +  T)  ©  s(t  +  T)  -  s(t  -  A  +  T)©s(t  +  T)  -DA(e) 

where  the  symbol  @  represents  modulo  2  addition.  Multiplication  of  ±1  terms  is 
equivalent  to  modulo  2  addition. 

In  order  to  determine  the  effect  of  the  self-noise  on  the  system  performance,  it  is 

necessary  to  determine  the  power  spectrum  of  n  (t  ,e  )  .  This  computation  is 

s 


The  symbol  E(  )  represents  the  ensemble  average  operator.  The  sequence  s(t) 

^is  assumed  to  have  an  equidistributed  random  time  origin. 

S.  W.  Golomb,  "Sequences  with  randomness  properties,"  Glenn  Martin  Co.  ,  Baltimore, 
gMd.  ,  Final  Report,  Contract  SC-54-36611,  June,  1955. 

B.  Elspas,  "A  radar  system,  based  on  statistical  estimation  and  resolution  considera¬ 
tions,”  Stanford  Electronics  Lab.,  Stanford  University,  Stanford,  Calif.,  Report 
#361-1,  August,  1955,  Appendix  D. 
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performed  for  values  of  delay  error  e  =  0  ,  ±mA  .  It  has  been  known  for  some  time 
that  maximal -length  linear  binary  sequences  possess  the  "cycle  and  add"  property,9, 10 

i.e.,  s  (t  +  iA)  ©  s(t  +  jA)  =  s  (t  +  rA)  .  This  relationship  can  be  used  to 
simplify  the  form  of  (4) : 


V»-£  =  mA>  s(t  +  T  +  jA)  -  a[(t  +  T  +  (j  +  1, a]  ,  for  m  =  0,  M,  elc. 

=  ±8(t  +  TtnA>  .  for  m  .  ±i.  Mil,  etc. 

s(t+T+rA)-s(t+f+qA)  ,  otherwise  (5) 

where  j  ,  n  ,  r  ,  q  are  integers  which  are  dependent  on  the  delay  error  e  =  mA  ; 
for  no  value  of  m  is  r  -  q  .  Hence,  it  is  seen  that  for  e  --  ±A  the  self-noise  power 
spectrum  is  the  same  as  the  signal  power  spectrum,  Gg(f)  .  The  signal  power  spec¬ 
tral  density  has  components  at  frequencies  which  are  integer  multiples  of  1/M  A  ,  i.  e.  , 


Vf>  =  MA  °(f> 


<5(f  t-  k/MA) 


l>  —  -oo 


where  6(f  fQ)  is  the  Dirac  delta  function  defined  by  the  operator  equation, 

J  H  ( f )  6  ( f  -  fQ )  df  =  H  ( fQ )  .  '  The  quantify  G  ( f )  ,  which  can  be  considered  as  the 
"envelope"  of  the  signal  spectrum,  is  plotted  in  Fig.  3a. 


i7 

1W-  Gol°mb,  "Sequences  with  the  cycle  and  add  property,"  Jet  Propulsion  Labor- 
DecTmber,  1957a  °f  Technol°^’  Pasadena«  Calif.  .  Section  Report  #8-573, 

10?,CD  Tits  worth,  L.  R.  Welch,  "Modulation  by  random  and  pseudo-random  sequences 

Jet  Propulsion  Laboratory,  California  Institute  of  Technology ,  Pasadena  California 
Progress  Report  #20-387 ,  June,  1959,  Sec.  3.  Pasadena,  California, 

In  other  words,  if  S  is  the  state  vector  of  the  FSR  at  the  time  t  +  iA  ,  and  T  is 
the  transformation  matrix  of  the  FSR,  then  (I  +  TJ  ~  i )  S  -  Tr  ~  i  S  where  I  is  the 

isUc  polynomial  o^T  "te  Tl"  -  \l\  '  V*"  r  f°,Und  Charactei- 

CjJey. -Hamilton  theorem  which  shows  that  0(T)  =  o”,  and  the  periodMtyrelation 
t  I  ,  to  obtain  the  relationship  between  these  integers. 
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Fig.  3  Envelopes"  of  power  spectral  densities,  (a)  binary  sequence  spectral 
density,  (b)  self-noise  spectral  density  for  e  =  A 
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G(f)  *  A 


For  other  values  of  delay  error,  the  self-noise  is  expressed  as  the  difference  between 
two  versions  of  the  same  sequence  having  different  time  origins,  and  the  self-noise 
spectrum  is 


Gns(f,e  =  mA)  =  2Gg(f)  [l  -  cos  27rfA]  •  ,  if  m  =  0  ,  M  ,  etc. 

r  (6) 

=  2Gg(f)  [1  -  cos  2rr  ( r  -  q)fA.]  ,  if  1  <  |m|  <  M  -  1  ,  etc. 

The  self-noise  spectrum  for  e  =  0  is  depicted  in  Fig.  3b.  Notice  that  it  is  a  general 
characteristic  of  the  spectra  in  (6)  that  they  have  a  null  at  the  origin,  and  as  a  result 
the  self-noise  is  more  easily  removed  by  the  loop  filter. 

m 

The  noise  input  to  the  discriminator  is  assumed  to  be  white  Gaussian  noise.  Because 
of  this  spectral  density,  it  produces  a  noise  component  in  the  output  of  the  cross- 
correlation  network  which  is  also  white 


k<5s(t  +  T)n(t)  =  kn  (t) 


(7) 


The  spectral  density  of  nn(t)  is  Gnn(f)  -  Pd  Nq  ,  where  Pd  is  the  average  power 
of  <5s(t  +  T)  and  Nq  is  the. input  noise  spectral  density.  For  values  of  T  which 
are  fixed  or  vary  little  in  a  time  interval  A  ,  the  average  power  in  6s  (t  +  T)  = 
s(t  +  T  +  A)  '-  s  ( t  +  T  -  A)  is  Pd  =  2.  Hence,  G^f)  =■  2Nq  watt-sec.  , 


Thus,  the  output  of  the  cross -correlator  network. can  be  written,  using  (1)  ,  (2)  ,  and 
(7)  as 


DA(f)  +  VVe)  +  nn(  t  >/ 


(8) 
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By  referring  to  Figs.  1  and  (8),  it  can  be  seen  that  the  system  equation  (in  operator 
notation)  is 


PT  *  F(p/po,  [da«t)  *  +  nn(t>//irj  (1 

where  gf  is  the  loop-filter  gain-constant,  pQ  is  the  loop-filter  frequency-constant, 
and  gc  is  the  gain  of  the  voltage-controlled  clock.  Define  gQ  as  the  dc  loop  gain, 
i.e.  ,  a  steady  delay  error  of  e  sec.  (\e\  <  A)  produces  a  gQc/A  cps  change  in' 
clock  frequency  or  a  gQe  sec.  delay  change/sec.  The  dc  loop  gain  is  given  by 
go  =  (M  +  !)/M  •  The  system  equation,  (9),  can  then  be  written  as 


pT  =  gQ  A  F  ( p/p 


■o>  (btt)  [d, 


'a(€)  +  V(t,£)  +  nn(t)/yp; 


(10) 


VI'  NOISE  PERFORMANCE 

In  this  section  the  effects  of  receiver  noise  and  self-noise  on  system  accuracy  are 

determined.  The  discriminator  is  assumed  to  be  in  the  locked-on  state,  i.  e.  , 

M  <  A  ,  so  that  the  system  is  operating  in  the  region  where  D  (e)  =  f(M  +  1)/Ml 
<fe/A).  ;  •  .  A  J 

Define  the  normalized  loop  gain  g  =  go/po  .  The  iinearized  system  equation  can  then 
be  obtained,  from  (10)  as 


■T/A  = 


F(P/P0) 

g'lp7ST 


£  +  Vl)/  \ZE 


A  +  (M  +  1  )/M 
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Since  e  =  T  -  T  ,  (10)  can  be  written 


a  /  n  (t  ,e)  +  n  (t )/./p~ 

T/A  -  H,p/po)  T/A  +  —  8 


where  H(p/pQ)  is  the  linearized  closed-loop  transfer  function 


.  a  SF(p/po) 

H(p/p0)  *  57p0+gF(p/po) 

The  linear  feedback  network  depicted  in  Fig.  4a  can  be  shown  to  have  the  same  system 
equation  as  (11),  and  therefore  it  provides  performance  equivalent  to  that  of  Fig.  1 

T/n  nD(t  ><0  +  n 

when  the  equivalent  input  —  ^ 1  ‘  +  — _ 11  . v  b  ,  ,  . 

.  4  cuemuipui  A  +  (M  +  1)/M  .  '  is  applied,  and  |e|  .<  A. 

#  . 

Consider  that  the  closed-loop  transfer  function  H(p/p  j  has  the  form 


H(p/p  )  = 


_  1 ' +  /2"  p/p0 

1  +  /T  p/p0  +•  (p/p0)2 


This  transfer  function  has  been  shown  optimum  for  ramp  inputs  of  delay  in  the  presence 
of  white  noise  in  that  it  minimizes  the  total  squared  transient  error  plus  the  mean 
squared  error  caused  by  interfering  noise. 12  The'  relationship  between  pQ  and  the 

transient  response  is  discussed  further  in  the  next  section!  The  loop-filter  shown  in 
Fig.  4b.  •  •  . 


npo?56f’  E'  R6Ch^’  "Design  arid  Performance  of  phase-locked  circuits  capable  of 

vol  'pt.imum  P®r^ormance  ovei'  a  wide  range  of  input  signal  levels,"  IRE  Trans. , 
vol.  IT-1,  pp.  66-72,  March,  1955. 
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F(P/P0) 


i  +  /sT  p/po 
1_+  &p/p0 


(12) 


can  be  used  to  approximate  this  closed-loop  response  as  closely  as  desired  by  making 
the  normalized  gain  constant  g  large. 


'Ilie  mean  square  delay  error  caused  by  the  receiver  noise  can  be  evaluated  from  (7) 
and  (11)  to  be 


P 

8 


G  ( f ) |  H  ( j  w/p  )|-2 
un  0 


df  *  2. 12A2 


-oO 


(13) 


For  example,  if  the  receiver  noise  temperature  is  900° K,'  pQ  =  l  rad/sec.  ,  M  »  1  ; 

A  =  10  6  sec.  and  Pg  =  10-15  watt,  then  Nq  =  6.  2  x  10"21  watt-sec. ,  and  the  rms 

delay,  error  caused  by  noise  is  an  =5.3  x  10“9  sec  ’ 

“n  ’ 

The  self-noise  spectrum  has  been  shown  to  depend  upon  the  delay. error  e  .  Hence, 
the  mean  square  value  of  delay  error  caused  by  the  self-noise  also  depends  on  e,  An 
upper  bound  on  the  mean  square  delay  error  can  be  found  by  assuming  the  most  pessi-. 
mistic  spectral  density  function,  i.  e.  ,.the  self-noise  spectrum  obtained  for  e  =  -±A". 
The  mean  square  self-noise  delay  error  is  then 


®  A 


/  S 


,(f)|H(jco/po)|2df  =  ^  S  ^J^/MAp0}  /sinjVM\2 

, l+(2^/MApQ)4  (  ,7rr/M  ‘ 


Under  conditions  where  the  number  of  digits  in  a  period  M  is  large  and  the  loop  filter 
constant  pQ  is  small  compared  to  the  signal  bandwidth  (~  1/2 A)  ,  this  expression 
simplifies  to 


(14) 
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For  the  parameters  just  considered,  p  =  1  rad/sec. ,  A  =  10  6  sec. ,  (14)  gives  an 

^  - 1 0 

rms  self-noise  delay  error  of  an  =  7.07  x  10  sec. 

"s 

When  the  delay  error  fluctuations  become  too  large,  the  discriminator  encounters  a 
threshold  effect.  This  effect  is  caused  by  the  fact  that  large  delay  errors  force  the 
discriminator  to  operate,  at  least  temporarily,  in  a  region  of  negative  loop  gain,  since 
D  (e  )  has  a  negative  slope  for  A  <  |  ej  <  2A  . 

Experimental  measurements  made  thus  far  indicate  that  the  system  threshold  occurs 
when  the  total  rms  delay  error  is  about  =  0.30  A  .  For  delay  error  having  a  < 
0.30  A  ,  experiments  showed  a  negligible  probability  of  losing  the  locked-on  state  in 
the  absence  of  transient  errors.  For  delay  error  fluctuations  which  have  a  Gaussian 
amplitude  statistic  and  o-^  =  0.  30  A  ,  the  probability  that  |  e|  >  A  in  the  linearized 
model  is  only  0.  . 00087.  This  statement  assumes  that  transient  errors  caused  by  fluc¬ 
tuations  in  T  are  small  compared  to  A  .  Assuming  that  the  self-noise  error  is  neg¬ 
ligible,  the  system  threshold  can  be  computed  from  (13)  and  is  given  approximately  by 


V  TRANSIENT  PERFORMANCE 

In  this  section  two  problems  concerning  the  transient  performance  of  the  discriminator 
in  the  absence  of  noise  are  investigated : 

•1.  How  rapidly  can  a  given  region  be  searched  and  the  target  acquired  for  a 

given  closed-loop  noise  bandwidth,  B  =  1.06  p  cps? 

n  Fo  F 

2.  Once  the  target  has  been  acquired,  what  is  the  maximum  change  in  velocity 
that  can  be  tolerated  without  losing  the  locked-on  condition? 
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In  order  to  solve  these  problems  most  readily,  we  rewrite  (10),  neglecting  the  noise 
effects,  and  making  use  of  the  following  time  normalizations  x  =  e/A  ,  y  =  T/A 
y  -  x  =  T/A  ,  r  =  pQt  ,  s  =  p/po  =  d/d t  ,  D ( x)  =  DA(e) 

s(y  -  x)  =  g  F  ( s  )  D  ( x )  (1 

The  integral  plus  proportional  control  loop-filter  of  (12)  now  becomes 

'  gF(s) 

Thus  we  obtain  the  operator  equation 

(l/g  +  s)s(y  -  x)  =  (1  +  yf%  s)D(x) 

This  system  equation  can  be  rewritten  in  time-derivative  notation 

y/g  +  y  =  x/g  +  X  +  D(x)  +  /2  D'(x)x  ^ 

where  x  =  dx/dr  ,  D'(x)  =  dD/dx,  etc. 

The  phase-plane  method  of  solving  this  second-order  differential  equation  is  to  com¬ 
pute  the  trajectories  in  x  ,  x  space  which  describe  the  solution  of  this  equation  for 
the  desired  sets  of  initial  conditions.  Define  a  new  variable  y  =  x/x  =  dx/dx 
which  from  (17)  is  given  by 

dx  '  r  ■  •  -o 

d?  -  nx,_x,  y  ,y] 


-  _  D(k)  +  \[2  D1  ( x)  +  l/gjx  -  y/g  -  y 


(18) 


_  1+^28 
1/g  +  S 
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Computer  solutions  to  these  trajectories  can  be  obtained  by  approximating  the  differ¬ 
ential  equation  (18)  with  the  difference  equation 


where  xq  ,  xq  are  the  initial  values  of  x  ,  x  ,  and  <5^  is  the  j-th  increment  in  x. 
Satisfactory  solutions  can  be  obtained  by  letting  the  computer  function  as  an  adaptive 

device  so  that  the  size  of  the  interval  is  made  variable.  In  this  particular  sequential 

13  14  ~ 

computation  the  interval  size  is  taken  as1 


6 

1  +  I  ,xj^ 


The  value  .6  to  be  used  is  0.02. 


■^R.  Bellman,  Adpative  Control  Processes:  A  Guided  Tour,  Princeton  University 
Press,  Princeton,  New  Jersey,  1961,  Chapter  4. 

Notice  that  the  distance  moved  in  the  .  x  ,  x  plane  in  one  increment  is 


Ar 

n 


-  x  ) 2  +  ( x  -  x  ) 2 
n  n  +  1  n; 


1  +  y2- 
'n 


1/2 


lI  +  2  y  +  y 

'  1  V  'iii 


Hence,  the  distance  moved. is  bounded  by  6// 2  <  r  <  <5  .  Thus,  the  maximum 
distance  moved  is  6  ,  yet  needlessly  small  increments  are  not  used  for  small  or 
moderate  [y|  . 
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Consider  the  search  and  acquisition  problem  where  the  normalized  loop  gain  g  =  »  . 
Assume  that  the  target  search  velocity  is  a  constant,  y(t)  =  y  ,  y  =  0.  The  variable 
y  now  becomes 


y(x,i)  ..muJLmxix 

x 

The  acquisition  trajectories  for  these  conditions  are  plotted15  in  Fig.  5.  If  the  delay 
error  is  decreasing  from  the  left  the  system  does  not  respond  until  x  =  -2.  As  can 
be  seen,  if  the  search  velocity  is  |y|  <  2.  2  ,  the  system  locks  on  and  the  state  vari¬ 
able  converges  to  the  origin.  As  an  example,  consider  A  =  lo-6  sec,  p  =  io  rps. 
Then  we  have  y  =  (dT/dt)/pQA  =  105  dT/dt  .  For  electromagnetic  propagation 
c  =  3  x  10  km/sec.  and  the  maximum  search  velocity  (two-way  propagation  time) 
permitted  is  v  =  (1/2)  2:2cpQA  =  3.  3  km/sec.  From  (15),  it  can  be  seen  that  the 
threshold  value  of  signal  power  is  Pg  =  6.36  Nqv/cA  =  70.  6  Nq  watt  in  this  example 
(M  »  1)  . 

Notice,  however,  that  even  if  the  system  fails  to  lock-on, in  this  pkrticular  interval  of 
x  ,  the  velocity  error  is  less  at  the  end  of  the  transient  than  at  the  beginning.  Since 
D  (x ).  is  periodic  every  M  ,  the  system  always  locks  on  eventually  because  with 
*  ■“  °°  there  is  no  decay  in  x  outside  of  the  intervals  |x|  <2.  It  may ,  however , 
pass  through  stable  regions  of  x  many  times  before  locking  on.  In  practice,  this 
behavior  might  rely  on  unrealistic  storage  times  in  the  filter,  as  is  implicit  in  the 
assumption  that  g  =  »  .  Furthermore,  the  time  required  for  lock-on  might  be  intol¬ 
erable  for  search  velocities  |yj  >'2.2. 


The  trajectoiy  computations  were  performed  by  Lt.  C.  S.  Mulloy  on  a  CDC  1604 
computer  and  are  presented  as  part  of  his  report  entitled  "Digital  analysis  of  the 
delay-lock  discriminator,"  U.  S.  Naval  Postgraduate  School,  Monterey,  California, 
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Fig.  5  Acquisition  trajectories  for  loop  gain  g  =  «,  plotted  for  various 
values  of  y  ,  the  normalized  search  velocity 


18 


LOCKHEED  MISSILES  &  SPACE  COMPANY 


The  effect  of  target  velocity  transients  can  also  be  ascertained  from  Fig.  5.  If  the 
system  is  initially  locked  on,  x  =  0  ,  x  =  0  and  the  target  suddenly  changes  velocity 
to  y  ,  the  system  response  is  described  by  that  portion  of  the  trajectories  which  be¬ 
gins  at  x  =  0  ,  x  =  y  .  As  can  be  seen,  normalized  velocity  transients  of  y  =  3.38 
can  be  tolerated  without  losing  the  locked-on  state. 

Acquisition  trajectories  have  also  been  obtained  for  finite  loop  gain  systems  with 
g  =  10  .  The  equation  for  y  in  this  example  is 

y(x,x,y)  =  - 

x 

Because  of  the  finite  loop  gain,  the'  maximum  steady-state  clock -frequency  change  is 
g  =  10 pQ  .  Thus,  for  this  reason  alone,  the  system  would  never  lock  on  for  ]y|  >  10. 

The  trajectories  for  these  conditions  are  shown  in  Fig.  6.  As  can  be  seen,  the  fact 
that  the  loop  gain  is  reduced  from  °°  to  10  has  very  little  effect  on  the  maximum  tol¬ 
erable  search  velocity;  [  y  [  £  2.2  is  still  permitted,  Notice  however ,  that  the  curves 

which  lock  on  converge  to  x  =  O.ly  or  e  =  xA  =  0.1  Ay  rather  than  to  x  =  0. . 

~6 

For  A  =  10  sec.  and  y  =  1  ,  the  steady -state  lock-on  point  is  e  =  0. 1  micro¬ 
seconds.  In  practice,  by  measuring  y  one  can  correct  for  this  steady-state  bias 
error. 

Although  it  is  not  shown  on  the  trajectory  curves,  the  value  of  x  for  the  curves  which 
do  not  lock  on  decays  with  a  time  constant  g/pQ  =  10/pQ  sec.  Thus,  if  the  period  of 
the  sequence  is  large  so  that  M  »  1  ,  the  value  of  x  will  have  decayed  almost  to  y 
when  x  =  M  -  2, 

As  well  as  knowing  whether  or  not  the  system  will  lock  on  for  a  given  search  velocity, 
it  is  also  important  to  know  how  long  the  transient  lasts.  The  time  response  of  the 
system  can  be  obtained  from  computer  solutions  to  the  difference  equation 


jy/2  D1 
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Fig.  6  Acquisition  trajectories  for  loop  gain  g  =  10  ,  plotted  for 
various  values  of  y  ,  the  normalized  search  velocity 


20 


LOCKHEED  MISSILES  8c  SPACE  COMPANY 


6-90-62-85 


t  +1  ~  t  =  t 
n  +  1  n 


(*» + i  5i)  -  i(xo * 2 6 i)  ®  — - 

\  i-o  /  \  i  =  0  /  (1  + 


I  y  | ) x(x  ) 
'  n ' '  '  n ; 


The  results  of  this  computation  are  shown  in  Figs.  7a,  b  as  plots  of  x  and  x  for 
g  =  00  ;  the  maximum  tolerable  search  velocity  is  y  =  2.  2  .  The  time  required  for 
the  transient  to  subside  within  J  xj  <0.1  is  j  =  5.  6  .  For  p  =  10  rad/sec.  this 
lock-on  time  is  about  0.  56  sec. 


VI  EFFECT  OF  AMPLITUDE- LIMITING  THE  RECEIVED  DATA 

Considerable  simplification  in  the  circuitry  is  possible  if  the  received  data  are 
amplitude-limited  before  entering  the  cross -correlation  network.  The  input  to  the 
cross -correlation  network  then  has  the  binary  form  u  ( t )  =  A  sgn  [^PgS  (t)  +  n{t)J 
where  A  is  the  limiter  output  amplitude  and  sgn  x  is  the  sign  function 

sgn  x  1  ,  if  x  >  0 
=  -1  ,  if  x  <  .0 

Thus,  the  waveforms  fed  into  the  discriminator  have  ±A  amplitudes,  and,  as  a  result, 
the  multiplication  circuits  of  Fig.  1  can  be  replaced  by  modulo  2  adders,  as  shown  in 
Fig.  8.  It  is  apparent  that  this  change  is  inconsequential  for  a  noise-free  input.  The 
purpose  of  this  section  is  to  determine  how  much  this  limiting  action  influences  the 
noise  performance  when  the  input  signal -to-noise  ratio  (SNR)  is  small. 

The  first  step  in  the  solution  of  this  problem  is  to  compute  the  cross-correlation  be¬ 
tween  the  output  of  the  limiter  u(t  +  p)  and  the  signal  component  s  (t)  .  It  is  this 
cross -correlation  which  determines  the  useful  signal  component  in  the  correlator  net¬ 
work  output  through  its  relationship  to  the  discriminator  characteristic.  Consider  that 
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DIFFERENTIATED  DELAY  ERROR  i 


or  g  =  oo  anr1 

~~l  delay  erro 


Fig.  7  Acquisition  transient  f< 
y  =  2.2  ,  (a)  normaliz 
delay  error  response 
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bandlimited  white  Gaussian  noise  enters  the  limiter.  The  noise  is  low-pass  having  a 
maximum  frequency  B  and  average  power  Pn  .  It  is  shown  in  the  appendix  that  the 
cross-correlation  is  related  to  the  signal  autocorrelation  function  by  the  equation 

E  [u(t  +  n)a  ( t )]  =  ARug(p)  =  ARs(p)  erf  yp^/2Pn  (2C 

where  erf  x  is  the  error  function.  For  small  input  SNR,  (20)  becomes 

n 

It  is  interesting  to  compare  this  cross -correlation  function  with  that  which  is  obtained 
without  the  limiter  when  the  average  total  input  power  is  held  constant  at  the  same 
level  A  .  This  type  of  operation  is  obtained  if  the  discriminator  is  preceded  by  an 
"ideal"  AGC  amplifier  which  has  constant  average  output  power  A2  .  The  cross¬ 
correlation  for  this  situation  is 


ARB(p) 


By  comparing  this  expression  with  (20),  it  is  apparent  that  at  small  values  of  input 
SNR^  the  cross -correlation  with  amplitude-limited  inputs  is  smaller  by  a  factor  of 
j/2A\_Since  both  expressions  are  proportional  to  the  square  root  of  the  input  SNR 
for  small  SNR,  the  loop  gain  of  the  discriminator,  gQ  ,  decreases  in  the  same  propor¬ 
tion.  More  generally,  the  loop  gain  with  limiting,  gQL  ,  varies  with  SNR  as 

goL  =  go  erf 

where  gQ  is  the  loop  gain  for  noise-free  inputs.  This  expression  is  plotted  in  Fig.  9. 
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It  is  next  necessary  to  determine  the  spectral  density  of  the  noise  output  of  the  cross - 

correlator  network.  It  is  known  that  at  low  SNR  the  autocorrelation  function  of  the 

1 6 

noise  in  the  limiter  output  is 


(tO 


?A2 

7T 


sin 


R  (ix) 
n'r; 

P 


n 


Hence,  the  noise  spectrum  in  the  limiter  output  is 


.  17 


G  If)  = 
no'  ’ 


2  A 
IT  P 


2  r 


n 


G  (  f )  +  0. 167  qG  ( f  )/P2  +  0.  075  gG  ( f  )/P4  +  . . . 
'  3  n'  7  n  5  n'  7  n 


Assume  that  the  noise  bandwidth  is  large  compared  to  the  signal  bandwidth.  Then  the 
noise  effects  in  the  cross -correlator  output  are  dependent  mainly  upon  Gno(0)  ,  and 
this  spectral  density  has  the  value 


G„„<°>  *  2  /  f  [lVM)/Pn]d’'  =  F  a2  /  8in‘‘ 

0  0 

(21) 

oo  f 

_  2  A2  f  sin  x  X  /sin  x\ 

v  tB  J  |  x  6  \  x  / 

0  L 


The  degradation  in  system  performance  at  low  SNR  caused  by  the  limiter  action  can 
now  be  determined  by  taking  the  ratio  of  the  square  of  the  signal  component  amplitude 


■*-6R.  Price,  "A  useful  theorem  for  nonlinear  devices  having  Gaussian  inputs,"  IRE 
Trans.  ,  vol.  IT-4,  pp.  69-72,  June,  1958. 

■^The  notation  3Gn(f)  represents  the  convolution  in  the  frequency  domain 
Gn(f)*Gn(f)*Gn<f)  . 
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at  the  cross-correlator  output  to  the  noise  spectral  density  Gfi0( 0 )  .  This  ratio  is 

then  compared  with  the  corresponding  ratio  obtained  with  an  ideal  AGC  amplifier  having 

2 

average  output  power  A  .  For  operation  with  the  limiter,  this  ratio  can  be  obtained 
from  (20),  (21)  as 


Bu>)/Rs%>  PS  /  1  \ 

Gno(0)  Pn/2B  \l.l/  <22) 

Operation  with  the  ideal  AGC  yields  a  ratio  which  is  larger  by  only  a  factor  of  1. 1  or 
0.  4  db.  Hence  the  use  of  an  amplitude  limiter  is  attractive,  since  it  allows  the  use  of 
binary  logic  for  the  cross -correlation  circuitry  at  only  a  small  expense  in  the  theoret¬ 
ical  noise  performance. 


VII  DISCUSSION 

The  transient  and  noise  analyses  show  that  the  delay-lock  tracking  system  possesses 
several  features  which  are  desirable  for  long-range  tracking  operations.  For  example, 
the  system  threshold  occurs  at  a  SNR  which  permits  operation  at  very  low  power  levels 
even  for  rather  large  target  velocities.  At  the  same  time,  the  signal  can  be  constructed 
to  have  periodicity  long  enough  to  resolve  any  range  ambiguities  which  might  otherwise 
occur  for  large  target  distances. 

Solution  of  the  phase-plane  trajectories  has  demonstrated  that  the  system  is  capable  of 
reasonably  large  search  velocities,  relative  to  the  loop-filter  time-constant,  which 
result  in  target  acquisition.  It  should  be  pointed  out  here  that  by  making  the  loop-filter 
adaptive,  certain  advantages  in  noise  performance  can  result.  One  can  sense,  by 
means  of  an  external  correlator,  whether  or  not  the  system  is  locked  on.  When  the 
loop  is  unlocked  and  the  discriminator  is  in  its  search  mode,  the  loop-filter  can  be 
designed  to  permit  target  acquisition  with  the  desired  search  velocity.  Once  the  target 


27 


LOCKHEED  MISSILES  &  SPACE  COMPANY 


6-90-62-85 


has  been  acquired,  the  closed-loop  bandwidth  can  be  decreased,  and  the  noise  per¬ 
formance  of  the  loop  can  thereby  be  improved.  More  extensive  discussion  of  this 
potential  will  be  presented  in  a  later  paper. 

Finally,  it  has  been  shown  that  amplitude-limiting  the  received  data  permitB  the  use 
of  binary  logic  in  place  of  analog  multipliers  at  only  a  small  expense  in  the  theoretical 
noise  performance.  In  practice,  deficiencies  in  performance  of  the  analog  circuitry 
would  often  more  than  make  up  for  this  theoretical  difference.  With  either  analog  or 
binary  inputs  to  the  discriminator,  if  the  average  input  power  is  fixed,  the  loop  gain 
decreases  in  proportion  to  the  square  root  of  the  SNR,  and  the  closed-loop  bandwidth 
and  maximum  search  velocity  are  diminished  accordingly.  Herice,  the  loop  gain 
should  be  made  large  enough  to  satisfy  the  search  velocity  requirements  at  the  lowest 
SNR  likely  to  be  encountered. 
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APPENDIX 


The  purpose  of  this  appendix  is  to  compute  the  cross -correlation  between  the  limiter 
output  u(  t  +  p)  and  the  reference  signal  component  s  ( t )  .  This  cross -correlation 
can  be  written 


E  z(t)  =  u-(t  =  r^)  =  Pr  [z  ( t )  =  1]-  Pr[z(t)  =  -1] 


where  u(t)  =  A  sgn  j^v/p^s(t)  +n(t)j  . 
s(t)  =  sQ  ,  s(t  +n)  =  s^  ,  n(t  +  p)  = 

1 )  =  p  =  ( M  +  1  )/2M  ,  and  Pr  { s  =  -1 )  = 
can  then  be  written 


In  order  to  simplify  the  notation,  define 
.  Define  the  signal  probabilities  Pr(s  = 
=  q  =  1  -  p  .  The  probability  that  z  =  1 


Pr  ( z  =  1 ) 


=  PPr<nM  >  '  =  !)  +qPr(n^  <  -  v^8/i|so  =  -1)  (23) 

=  Pr  (n^  >  -  /F)[pPr(s^  =  l|so  =  1)  +  qPr(s^  =  -i\bq  =  -l)j 

+  Pr(n^  >  /P^jpPr^  =  -l|sQ  =  1)  +  qPr(s^  =  1[sq  -  -l)j 


where  use  has  been  made  of  the  fact  that  the  probability  density  of  n  is  symmetric 
about  the  origin.  The  probability  that  z  =  1  can  be  obtained  in  a  similar  manner 

Pr(z  =  -1)  =  Pr  (n^  <  -^)  ^pPr(s^  =  1 1  sq  =  1)  +  qPr(s^  =  -l|  sq  -  -l)j 

(24) 

+  Pr(nM  <->/Pg)  [pPr(sM  =  -ll  sq  =  1)  +  qPr(s^  =  l|  sQ  =  -1) 
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The  autocorrelation  function  of  the  signal  is  expressed  by 


Rg(p)  =  E  [s(t)s(t  +  m)]  =  P  Pr(s^  =  l|so  =  1)  -  Pr(s^  =  -l|  sq  =  1) 


+  q  Pr ( =  -l|so  =  -1)  -  Pr (s^  =  l|sQ  =  -1)J 


Hence,  by  using  (23),  (24),  the  cross -correlation  function  R  (p)  can  be  written 

UB 


Rus(">  *  V*)  Pr<IV 


Under  the  assumption  that  the  noise  has  stationary  Gaussian  amplitude  statistics  and 
has  a  mean  square  value  Pfi  ,  we  obtain 


R  (p)  =  R  (M)  erf  /P  /2P 
us  s  v  s  n 
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